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ABSTRACT 

GSH 91.5 + 2 - 1 14 is a large HI shell located in the outer Galaxy at a kinematic distance of about 15 kpc. It was first identified in the 
Canadian Galactic Plane Survey (CGPS) by Pineault et al. (2002) as being possibly associated with objects possessing infrared colors 
which indicates strong stellar winds. The HI shell has no obvious continuum counterpart in the CGPS radio images at 408 and 1420 
MHz or in the IRAS images. We found no evidence for early-type massive stars, most likely as a result of the large extinction that 
is expected for this large distance. An analysis of the energetics and of the main physical parameters of the H I shell shows that this 
shell is likely the result of the combined action of the stellar winds and supernova explosions of many stars. We investigate whether a 
number of slightly extended regions characterized by a thermal radio continuum and located near the periphery of the H I shell could 
be the result of star formation triggered by the expanding shell. 
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1. Introduction 

A massive star possessing a strong stellar wind (SW) can inject 
as much energy into the interstellar medium (ISM) during its 
lifetime as it does during the final supernova explosion. If the 
stellar spatial velocity with respect to the ambient ISM is not 
too large (i.e. < 30 kms -1 ), the SW is expected to evacuate a 
large cavity around the star (called a stellar bubble) surrounded 
by a ring or shell of enhanced density. These structures have 
been observationally detected as optical and infrared (IR) nebu- 
; lae around Wolf-Rayet (WR) and Of stars (??). 

The situation with respect to radio observations, reviewed 
by ?, has evolved considerably with the advent of large Galactic 
neutral hydrogen (HI) surveys, namely the Canadian Galactic 
Plane Survey (CGPS; ?), the Southern GPS (SGPS; ??) and the 
VLA (VGPS; ?). These surveys at ~ 1 arcmin resolution have 
made it possible to identify and study the structure and dynamics 
of many neutral H I shells (????). Interestingly enough, the large 
majority of shells detected by their neutral hydrogen emission 
have low inferred expansion velocities, typically less than or on 
the order of 10 km s _I (?). 

In parallel to these developments on the observational scene, 
recent theoretical studies, building up on the initial work of ? and 
others, have considerably increased our understanding of the in- 
teraction of stellar winds with their surrounding ISM. The effects 
of the different evolutionary phases and of a large peculiar mo- 
tion of the star have been modeled in detail in a number of new 
studies (e.g., ????, and references therein). Concerning the often 
observed low expansion velocities of ~ 10 km s , ? have shown 
that velocity dispersion within the shell and the role of the local 
ISM background may significantly affect the appearance of an 
expanding HI shell in velocity space and thus its inferred pa- 
rameters (in particular, mass and expansion velocity). 



In the case of moderate-size shells, a puzzling aspect is the 
apparent lack of a radio continuum counterpart or of a candi- 
date progenitor star (????). This suggests that either one or sev- 
eral basic ingredients are missing in the predictions of the theory 
and/or that some detected H I shells are not real. 

Increasing the sample of well studied H I shells is a first natu- 
ral step in elucidating some of the current puzzles. In an attempt 
at diversifying the sample of known H I shells (consisting mostly 
of shells around objects known for their optically interesting fea- 
tures), ? initiated a project aimed at discovering new SW source 
candidates by using first IRAS colors to extract potential can- 
didates and then the CGPS database to look for a morphology 
indicative of a SW, i.e. shells, rings, bubbles, cavities, or voids. 
An obvious advantage of this procedure is that potential candi- 
dates suffer much less from the selection effects associated with 
optically chosen targets, for example, WR or Of stars, the distri- 
bution of which is severely biased by absorption. A positional 
coincidence between one or more candidates and a shell-like 
morphological structure is nevertheless not a proof of a phys- 
ical association, and a detailed analysis is required before any 
firm conclusion can be drawn. 

In this paper, we focus our attention on a very large (nearly 
1°5 diameter) and symmetrical H I shell centered at (I, b) = (91°5, 
+2°.Q) and observed at a velocit)0 v\ s[ » -114 kms -1 , which 
was identified by ?. A flat rotation curve model for the Galaxy 
gives for the shell a distance d = 15 kpc, a galactocentric radius 
R g = 17 kpc, a distance z ~ 525 pc from the Galactic plane, and 
a diameter D = 400 pc, placing it in the outermost part of the 
Galaxy. 

This shell thus offers the opportunity to explore the envi- 
ronment of remote regions of the Galaxy where many physi- 



1 All velocities are with respect to the local standard of rest (lsr) 
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Table 1. Observational parameters for the H I data 



Parameter Value 

Synthesized beam V.26 x 0198 

Observed rms noise (single channel)(K) 1.6 

Bandwidth (MHz) " 1.0 

Channel separation (km s -1 ) 0.824 

Velocity resolution (km s" 1 ) 1.32 

Velocity coverage (km s" 1 ) 224 



cal parameters such as metallicity, density, smaller or negligi- 
ble perturbations from spiral arms, greatly differ from those in 
nearby regions of the Galaxy. Despite this remoteness, there is 
evidence that star formation in the outer Galaxy may be com- 
mon, as shown by the discovery of a considerable number of 
embedded star clusters in molecular clouds up to galactocen- 
tric radii R g « 17kpc (??). A particularly interesting case 
is the suggestion by ? that star formation in Digel's Cloud 2 
(R g w 19 kpc) could have been triggered by the huge supernova 
remnant GSH 138-0 1-94 previously discoveredby ?. These huge 
shells imply that massive stars form in the outer Galaxy, empha- 
sizing the importance of studying these objects. At these large 
distances, optical obscuration is very severe, so that one has to 
resort to radio or infrared observations. 

The plan of the paper is as follows. In Section[2]we describe 
the observational data used, in Section [3] we briefly review the 
initial data used by ? and present detailed neutral hydrogen (H I) 
and continuum images of the object. The results are analyzed 
and discussed in Section [4] Section [5] is a summary of the main 
conclusions. 

2. Observational data 

Radio continuum data at 408 and 1420 MHz and 21 -cm spectral 
line data were obtained at DRAO as part of the CGPS survey 
(?). A detailed description of the data processing routines can 
be found in ?. At the position of GSH91.5+2-114, the contin- 
uum images have a resolution of 1 '08 x 0.'82 and 3'65 x 2f 82, 
and a measured noise of 0.082 K and 0.75 K, at 1420 and 408 
MHz respectively. Parameters relevant to the H I data are given 
in TableQ] High-resolution-processed (HIRES) (?) IRAS images 
produced at the Infrared Processing and Analysis Center (IPAC j3 
were also used. The images used are the result of 20 iterations 
of the algorithm, giving an approximate resolution ranging from 
about 0:5 to 2 '. At the position of GSH 91.5+2-1 14, the 60 mi- 
cron image has an approximate resolution of 1 '. 

3. Results 

3.1. The initial data 

Figure [TJ shows an H I image averaged between v = -1 10.3 and 
-1 16.9 km s _1 , showing the large H I shell discovered by ? with 
the four IRAS sources which they suggested might be physically 
associated. 

The interior of the shell seems to have been entirely cleared 
of neutral hydrogen. Though the shell is generally quite well de- 
fined, its northern part is essentially absent, or at least it does not 
show a clear outline. Indeed, the general morphology suggests 
that the shell is open to the north in a direction away from the 

2 The Infrared Processing and Analysis Center (IPAC) is funded by 
NASA as part of the Infrared Astronomical Satellite (IRAS) extended 
mission under contract to the Jet Propulsion Laboratory (JPL) 
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Galactic Longitude 

Fig. 1. HI emission distribution averaged between 
-110.3 and -116.9 kms -1 . The plus symbols indi- 
cate the position of the four DWCL source candidates, 
IRAS 21 147+5016, IRAS 21 106+5013, IRAS 21093+4959, 
and IRAS 21041+5052. The inset is a close-up view of the 
region near IRAS 21 147+5016. 

Galactic plane. Indicated on Fig.QJare two filaments (denoted A 
and B) that could be related to the shell, forming an incomplete 
northern border. A closer inspection of the figure also shows that 
the H I shell is far from homogeneous. Moreover, there are many 
quite well defined cavities visible that are projected onto the H I 
ring. 

The position of the four IRAS sources is indicated by plus 
symbols in Fig.QJ Two of the sources are seen projected well in- 
side the shell structure, while IRAS 21093+4959 is located be- 
tween filaments A and B. IRAS 21 147+5016 is seen projected 
onto a smaller H I cavity centered at (I, b) = (92°. 1, 0°9) and lo- 
cated on the larger shell structure. The inset in the bottom right 
corner of Fig. QJshows this region in more detail. Table|2]shows 
the IR colors of these four sources. 

A comparison of these colors with the values given by ? in 
his Table 4 Q allows us to conclude that the IRAS sources are 
probably dust shells related to WC8-9 stars (dubbed dusty late- 
type WR or DWCL by Cohen). In their preliminary analysis of 
the then-incomplete CGPS, ? had concluded that the surface dis- 
tribution of the stellar wind candidates chosen on the basis of 
their IR colors was on the average of 1 source per 10 square 
degree area. Note that there are no other sources with the same 
properties in the entire field shown in Fig. QJ The position of 
the four candidate sources inside or near the boundary of the H I 
shell prompted a more detailed analysis. 

3.2. Single-channel H I 21 -cm data 

Figure [2] shows the H I distribution in the LSR velocity range 
from -107.0 to -121.0 kms -1 . The original images were 

3 Cohen's (1995) definitions are of the form [i] - [J] = kjj + 
2,51og(F ; -/F i ) where the constants fey are 1.56 and 1.88 for the colors 
[12] - [25] and [25] - [60], respectively. 
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IRAS source 


/ 


b 


C n 


c 23 


F k 


(Jy) 


L ir (10 j L G ) 


21147+5016 


92° 05 


0? 97 


-0.34 


0.21 


9.5 


±0.2 


3.0 ± 1.2 


21106+5013 


91° 56 


1°41 


-0.32 


0.02 


11.2 


±0.2 


4.0 ± 1.6 


21093+4959 


91° 25 


1°39 


-0.34 


0.16 


9.2 


±0.2 


3.3 ± 1.3 


21041+5052 


91° 34 


2° 60 


-0.32 


-0.33 


27.1 


±0.2 


9.6 ±3.8 



DWCL refers to a dusty late-type Wolf-Rayet star. 

Cjj = log(Fj/Fj), where i and j run from 1 to 4 and correspond to 12, 25, 60 and 100 /jm, respectively and F, is the flux density (Jy) in band i. 
Ljj/L = 1.58 Fj r c , where the integrated flux in Jy is F„ = 1.3(Fi2 + Fis) + 0.7(F25 + F6o) + 0.2(F6o + -Fkx)) (?)■ Subscripts are wavelengths 
in microns . Luminosity values are for d = 15 kpc. 



smoothed to a 3-arcmin resolution to increase the signal-to-noise 
ratio. At v = -109.5 km s _I the shell, which we shall now refer 
to as GSH 91.5+2-1 14, is clearly observed and it is well defined 
down to v = -1 18.5 km s _1 . This structure is centered at (I, b) = 
(91°5, +2°) and has an angular size of about V.5. 

A limited region of bright centrally-located emission at - 
107.0 kms 1 could form part of the so-called receding cap, im- 
plying an expansion velocity of some 7 kms -1 . In an ideal situ- 
ation, one would then expect to find an approaching cap at a ve- 
locity of about -121 kms _I . No obvious excess HI emission is 
seen near this velocity. The lack of confusing emission at this ve- 
locity suggests that either expansion on the near side took place 
in an extremely low-density medium, making its detection be- 
low the sensitivity of the telescope, or that the shell was very 
incomplete and nearly absent on the near side. 

3.3. Radio continuum and infrared data 

? had already noted that no continuum (radio or IR) counterpart 
seemed to be present. Figure[3]shows the CGPS 1420 MHz radio 
continuum and 60 /mi HIRES images of the same field of view as 
the preceding figures. Discrete point sources have been removed 
from the radio continuum image and the shadings chosen so as to 
highlight the low-level emission centered on GSH91. 5+2-114. 
The relatively bright and large incomplete shell-like structure to 
the west of GSH 91.5+2-1 14, near (/, b) = (90°9, 1°6) is the H II 
complex BG 2107+49 discussed by ? and located at a kinemat- 
ical distance of 10 kpc. The circles delineate the approximate 
inner and outer boundaries of the H I shell, revealing that there 
is no obvious radio continuum or infrared emission related to 
GSH 91.5+2-1 14. 

The faint radio continuum emission seems to consist of fila- 
ments more or less aligned perpendicular to the Galactic plane, 
especially toward the north of the image, above the opening in 
the H I shell. The diffuse emission is more or less distributed ho- 
mogeneously over the northern part of the image, although there 
appears to be a slight excess of emission within the shell bound- 
ary, in particular in its northernmost part. In the south, the loca- 
tion where radio emission becomes fainter corresponds roughly 
with the H I shell boundary. 

Two interesting and slightly extended (6' - 8') structures 
(named G9 1.56+0. 97 and G92.24+1.57) are seen near the inner 
boundary of the HI shell. A close-up view of these structures 
is presented as insets at the bottom of Fig. [3] These objects are 
discussed in more detail in Section l4~4l 

We also examined the CGPS polarization images (percent- 
age polarization and position angle) for this field and did not find 
any evidence for highly structured emission associated with the 
shell or any of the slightly extended structures discussed above. 
This agrees with the findings of ?, who concluded that all polar- 



ization structures observable in the CGPS are generated and/or 
Faraday rotated closer than 2.5 kpc. 

4. Discussion 

Some of the questions we shall try to answer are: can the ob- 
served HI shell be a stellar wind shell formed by the winds 
from one or more of the four candidate DWCL sources? If 
not, what are the alternative shell-formation scenarios? Was 
the smaller HI cavity (the one embedded in the wall of the 
larger HI shell) caused by IRAS 21 147+5016 and, if so, could 
IRAS 21 147+5016 represent a case of triggered star formation? 
In this section we shall explore different possible formation sce- 
narios based on the determined size, mass, age, and energy in- 
volved. From here on we assume a distance of 15 +3 kpc for 
GSH 91.5+2-1 14, as suggested initially by ?. 

We note however that this distance is likely an upper limit. 
Indeed, using a new method (?) based on HI column densi- 
ties for the determination of distances within the disk of the 
Galaxy, ? obtained a new distance for the outer Galaxy HII 
region CTB 102. As this object is at I = 93° 115, b = 2°835, 
i.e. very nearly along the same line of sight as our shell, we 
can use ?'s Fig. 5 to extrapolate their velocity-distance rela- 
tion to obtain a distance estimate of approximately 1 1.5 kpc for 
GSH 91.5+2-1 14. This slightly lower distance still corresponds 
to an object well in the outer Galaxy at R g » 14.5 kpc and z ~ 
400 pc. Whenever possible, we shall show the distance depen- 
dence of simple parameters by defining dis - d/\5 kpc. 

4.1. Parameters of the H I shell 

Error estimates for all the parameters are based on our ability 
to estimate the spectral extent of the shell in velocity space, its 
spatial extent, and distance. 

The HI ring is clearly observed over some 14 kms -1 . 
Following the procedure described by ?, we derive a mass for 
the shell of M sh = (2.3 + 1.3) x 10 5 d 2 5 m Q . Note that this 
value agrees with that obtained using ? equation for the mass 
swept up by a shell, namely M s h(m ) = 8.5/?^ h (pc) which, tak- 
ing R sh = (200 + 30) dis pc, gives M sh = 3.4 X 10 s d\ 5 m Q . Using 
our determined value for the mass and considering a spherical 
cavity, we obtain for the initial (i.e. before the gas was swept up 
in the shell) neutral gas particle density n„ = (0.3+0.1) d^ cm -3 . 

Adopting an expansion velocity equal to half the velocity 
interval where the structure is observed, v exp = 7+2 kms -1 , 
we obtain for the kinetic energy stored in the expanding shell 
£kin = M sh v 2 exp /2 = (1.1 +0.6) x 10 50 ^ s erg. The kinematic age 
of the shell is given by f(Myr) = a/? s h(pc)/v eX p(kms _1 ), where 
a — 0.25 for a radiative SNR shell and a = 0.6 for a supershell. 
We obtain t = (7±2)di 5 Myr(a = 0.25) and t = (17+5) d X5 Myr 
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Fig. 2. CGPS gray-scale images showing the HI distribution in the LSR velocity range from -107.0 to -121.0 kms -1 . Velocity 
resolution is 1.3 km s -1 , and spatial resolution is 3'. The LSR central velocity of each image is indicated in the top left corner. The 
star symbols indicate the position of the four DWCL source candidates. For all images, darker shading indicates higher brightness. 

(a = 0.6). With these determined observational parameters as a association at all and that what we observe is simply a pure line- 
basis, we are now in a position to discuss the different formation of-sight coincidence, 
scenarios for GSH 91.5+2-1 14. 



4.2.1 . The action of the DWCL source candidates 



4.2. Formation scenarios for the H I shell 



The long kinematic timescale of 17 Myr, if GSH91.5+2-1 14 
has a SW origin, suggests that more than one generation of mas- 
In this section we assume for simplicity that the DWCL objects sive stars should be involved and/or that contributions from one 
are located at the same distance as GSH 9 1.5+2-1 14, namely 15 or more SN explosion are to be expected. Nevertheless, as a first 
kpc. Of course we cannot rule out that there may be no physical step, it is of interest to estimate the contribution that the winds 



S. Cichowolski & S. Pineault: A large HI shell in the outer part of the Galaxy 
VEL{LSR): - 1 1 4.4-1 km/e VEL(LSRJ: - 1 1 5.24 km/B VEL{LSR): -115.06 km/B 




J I L 



VEL(LSR): -1 16.63 km/s 



VEL(L5F!); -1 19.35 km/s 



J I I L 



* 



J I I L 



VEL(L£RJ: -117.71 tm/s 



VEL{LSR): -1 1S.&+ krn/s 




~i 1 1 r 



r 

r* * * 



J I L 



n 1 1 r 



* * 



J I L 



VE^LSR); -121. 01 km/s 




i 1 r 



* * 



j i i u 



"i 1 r 



+ * 



j i i i_ 



93° 92° 91° 90°- 

GALACTIC LONGITUDE 



Fig. 2. cont. 



of one or more of the DWCL candidates could make to the for- 
mation of this large H I shell. 

We first calculated the total integrated IR luminosity of each 
candidate. The derived parameters are shown in Table|2] If, fol- 
lowing ?, we assume the objects to be WC8-9 stars, then their 
total luminosity is about (?) L/L Q — 1.1 x 10 5 . The IR luminosi- 
ties thus represent only about from 3 to 10 % of the total stellar 
luminosity. 

Ignoring any intrinsic absorption by dust surrounding 
DWCL stars, we may ask whether a WC8-9 star would be eas- 
ily identified at a distance of 15 kpc along the line of sight to 
GSH 9 1 .5+2- 114. From the foreground H I column density and 



distance we can calculate typical values for the extinction and 
reddening of stars. The measured foreground H I column density 
is about 1.1 x 10 22 cm -2 , which corresponds to a reddening of 
E B - V = 2.3 mag (E B -v = AW4.8 x 10 21 cirr 2 ; ?). This gives a 
visual absorption A,, of 7.3 mag. Taking for the absolute magni- 
tude of WC8-9 stars (?) M v = -4.5, we obtain m v = 18.7. A star 
of this magnitude would clearly not stand out among field stars, 
although it would be detectable. In their study of GSH 90+03-99, 
? had estimated a visual extinction of 6.4 mag toward / ~ 90° and 
concluded that deep measurements were needed to detect early- 
type stars at a distance, for their object, of 13 kpc. 
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Fig. 3. 1420 MHz radio continuum and 60 yum images from the CGPS survey. Resolution for both is about V. Point-like sources 
have been removed from the radio continuum image. Circles indicate the approximate inner and outer boundary of the H I shell seen 
in Figure Q] The two insets at the bottom of the radio image show a close-up view of G9 1.56+0. 97 and G92.24+1.57 (point-like 
sources included, see text). 



There are reasons however to believe that the value of ab- 
sorption derived from the HI column density is a significant 
underestimate of the true value. As an alternative, we used the 
2MASS Point Source Catalogue (?) to produce a H - K s vs 
J - H diagram (Figure 0J. There are 1573 2MASS sources in 
a circular area centered at (I, b)-( 9V.1, 1°83) within a radius of 
400". These sources are indicated by green dots in Fig. [4] The 
positions of the dereddened early-type main sequence and gi- 
ant stars are indicated with blue and red solid lines, respectively. 
The blue and red dashed lines show the reddening curve for 09 
V and M0 III stars, respectively. Thus, normally reddened main 
sequence stars lie between the two dashed lines. Because absorp- 
tion is proportional to distance, we can infer that statistically the 
most distant sources have a visual absorption of about 13 mag 
along this direction in the Galaxy. Hence, as we are dealing with 
a structure located in the outer part of the Galaxy at about 15 
kpc, it is likely that the visual absorption is significantly higher 
than 7.3 mag. If absorption is as high as 13 mag, the expected ap- 
parent magnitude of WC8-9 stars could be as faint as m v = 24.4. 
The higher absorption derived from the color-color diagram sug- 
gests the presence of molecular gas along the line of sight in this 
direction. It also makes the optical identification of any massive 
star at that distance even more problematic. 

Note that lowering the distance to 11.5 kpc only makes the 
above magnitudes brighter by 0.6 mag. 

Next we estimate the energy injected by the stellar winds, 
E w = Mv\, t/2. Adopting (?) for a single WC8-9 a mass loss 
rate of M = lO" 5 m yr _1 , appropriate for a solar metallic - 
ity Z , and a wind velocity of v w = 1500 kms -1 , we obtain 
E w = 2.25 x 10 44 f(yr) erg. If the WC phase lasts ~ 10 5 yr, each 
star would impart E w = 2.25 x 10 49 erg to its local ISM during 
this phase. But only a fraction e of this energy gets transferred 
to the gas. According to evolutionary models of interstellar bub- 
bles, the expected energy conversion efficiency e = E^ n /E w is 



on the order of 0.2 or less (?). Observationally, the values of e 
derived from optical and radio observations can be as low as 0.02 
(?), indicating that there are cases where severe energetic losses 
may occur. 

Furthermore, given the location of GSH 9 1.5+2- 114 in the 
outer Galaxy, the effect of metallicity should be taken into ac- 
count. ? suggested that the mass-loss rate depends on metallicity 
Z as M oc Z'", with m = 0.5. The velocity of the wind is also 
lower at lower metallicity (?). ? concluded that the wind energy 
E w decreases by a factor of about 3 when Z decreases by a factor 
of 10 for stars with the same luminosity. Thus, for Z - 0.1 Z Q , 
each WC8-9 star would simply inject E w ~ 7.5 x 10 48 erg in 
the ISM during the WC phase. This implies that even with a 
relatively high efficiency e = 0.2, the injected kinetic energy is 
only £kin = eE w = 1.5 x 10 48 erg. Summarizing, neither for 
Z = Z nor for Z = 0.1 Z could the HI shell have been cre- 
ated only by the stellar winds of the DWCL sources in their WC 
phase. Even by considering a comparable contribution from the 
previous evolutionary phase of each star, the available kinetic en- 
ergy falls considerably short of the observed shell kinetic energy 
Skin = (1-1 ±0.6) xl0 5 X 5 erg. 

4.2.2. The action of unseen massive stars 

We have failed in finding any massive star in the area. However 
we have seen in Section l4.2.1l that distance and absorption con- 
spire to make even bright early-type stars inconspicuous. At a 
distance of 15 kpc and with an absorption of 7.3 mag, based 
on the HI column density, an 03V star (M v = -6.0) would 
have m v = 17.2, whereas a B1V star (M v = -3.2) would 
have m v = 20.0. For an absorption of 13 mag, derived from 
the 2MASS color-color diagram, the corresponding magnitudes 
would be 22.9 and 25.7, for an 03V and a Bl V star respectively. 
As mentioned above, a shorter distance of 1 1 .5 kpc leads to es- 
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H-Ks 



Fig. 4. (H-K,J-H) color-color diagram of the 2MASS sources 
found in a 400"circular area centered at (l,b) = (91°7, 03). 
The two solid curves represent the location of the main sequence 
(blue line) and the giant stars (red line) derived from ?. The par- 
allel dashed lines are reddening vectors with the crosses placed 
at intervals corresponding to five magnitudes of visual extinc- 
tion. We assumed the interstellar reddening law of ?, Aj/Ay = 
0.282 and A H /A v = 0.175. 

timated magnitudes brighter by only 0.6 mag. Although one or 
more early-type star cannot be easily detected at the distance of 
GSH 91.5+2-1 14, we cannot rule out the presence of these stars 
and their possible role in forming the HI shell. Furthermore, 
given that massive stars evolve in a short time, the originating 
stars might have vanished a long time ago and be not observable 
any more. 

If the origin of GSH91.5+2-114 is indeed the action of 
many massive stars, from the estimated value of the kinetic en- 
ergy of the shell (E kin = (l.l±0.6)xl0 50 erg, see Section l4~TT i we 
can estimate how many massive stars would have been needed to 
create it. Adopting mean stellar wind parameters for O-type stars 
(i.e. M = 2 x 10~ 6 m yr" 1 and v w = 2000 km s -1 , ?) and consid- 
ering that the main sequence phase lasts for at least 3 x 10 6 yr, the 
injected wind stellar energy is about 1.4 x 10 50 erg. Bearing in 
mind that the conversion efficiency is about or less than 0.2 and 
considering the possible effects of a lower metallicity environ- 
ment, the energy released by the stellar winds should be greater 
that 5.5 x 10 50 erg. Thus, at least four O-type stars were neces- 
sary to create GSH91. 5+2-114. In the case of a lower distance 
of 1 1 .5 kpc, the energy released by the stellar winds should be 
greater than 3.2 x 10 50 erg, implying that at least three O-type 
stars were required. 

4.2.3. A supernova explosion 

? proposed and interesting procedure for setting various con- 
straints on an SNR origin, which they applied to their H I shell 
and another one previously discovered by ?. Their procedure es- 
sentially rests on the expression giving the maximum observable 
radius of an SNR before it merges with the ISM (?): 

Emerge — 51 .3 Eq YIq Z ^ pC, 

where Eq is the explosion energy in units of 10 51 erg, no 
the mean ambient particle density in cm -3 , and Z the metal- 
licity normalized to the solar value. The input parameters to 



the above equation are the observed radius and particle den- 
sity, /? s h = 200c/i5pc and no = 0.3 d^ cm 4 . The results for 
GSH91.5+2-114 are shown as the last entry in Table[3] where 
the numerical values are for d = 15 kpc. The first two entries are 
identical to Table 2 of ? (with the addition of two parameters and 
a slight change in notation for added clarity). 

The column labelled R* mge (column 5) gives the merging ra- 
dius assuming a canonical value of Et = 1 for the explosion 
energy. As with the other two shells, the observed radius /? s h 
is about twice as large as the predicted merging radius R* melge , 
which means that the shell should have vanished well before 
reaching the observed radius. This implies that either Eq is larger 
and/or the object is closer. The column labelled Eg (column 
6) gives the energy required for R melge to be equal to the ob- 
served radius of the shell. The morphology of GSH 91.5+2-1 14 
is clearly too regular and well defined to represent a merging 
shell so that the value of E^ is a lower limit. The last column 
is the maximum distance GSH 9 1.5+2-1 14 would have to be in 
order to be the result of a single SN explosion with Eq = 1 ■ This 
distance corresponds to a systemic velocity of about -52 km s _1 
which, as is the case for the other two shells, is totally unrealistic, 
given the observed systemic velocity of -1 14 km s 4 . 

If we compare the three shells, we see that the first and last 
one are very similar in radius and mass, because they are vir- 
tually at the same distance and evolving in an ISM of com- 
parable particle density. The second one, GSH90+03-99, is 
somewhat closer and appears to be located in a significantly 
denser medium. Nevertheless, as discussed by ?, neither of the 
three shells could have been caused by a single SN explosion. 
Furthermore, as emphasized by ?, this argument applies as well 
to a single O- or B-type stellar wind bubble. Using 11.5 kpc as 
the distance would only change R„ er g e an d ^max by about 10% 
and E* by about 40%, thus not altering the general conclusion. 

4.3. Additional constraints from radio continuum emission 

From Figure [3] it is clear that there is no evidence of enhanced 
radio continuum emission associated with GSH91. 5+2-114. 
Although this absence is somewhat puzzling for H I shells as- 
sociated with known optical objects such as WR or Of stars, it is 
less so here. The kinematic age of the shell could be as high as 17 
Myr and massive O stars and early B stars could well have ended 
their lives. However, because of the large distance and large ex- 
tinction involved, we cannot completely rule out the presence of 
one or more of these massive stars. This statement can however 
be made more precise. 

Bellow we shall use the following parameters for the HI 
shell (note that 1 ' is equivalent to 4.5 pc at a distance of 15 kpc): 
inner and outer radius, = 160pc, R = 200 pc, shell thickness 
AR = 40 pc and initial particle density n = 0.3 cm 4 . These are 
average quantities because obviously the shell radii and thick- 
ness vary significantly with azimuthal angle. With these param- 
eters, we deduce that the particle density within the shell is now 
n s - 2.6 x n„ — 0.8 cm 4 . For simplicity, we shall neglect the 
presence of helium in the following simple estimates. 

4.3.1 . The simplest case: a homogeneous shell 

We assume that an ionizing star is located at the center of the H I 
shell and ask to what extent the inner side of the H I shell should 
be ionized, assuming that negligible material is present inside 
the shell. This is the classical Stromgren sphere problem where, 
instead of integrating from to some maximum radius R m , we 
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Table 3. Constraints on a single SNR origin 



Object 


«o 




Z 


merge 


e; 




Reference 




(cm- 3 ) 


(pc) 


(Zo) 


(pc) 


(10 51 erg) 


(kpc) 




GSH 138-01-94 


0.27 


180 


0.1 


93.3 


8.0 


5.9 










1.0 


83.0 


11.6 


4.9 




GSH 90+03-99 


2 


55 x 110 


0.1 


44.7 


6.3 


5.2 


2 








1.0 


39.8 


9.1 


4.3 




GSH91.5+2-114 


0.3 


200 


0.1 


90 


12.6 


6.9 


This paper 








1.0 


80 


18.2 


6.1 





^meigc (column 5): merging radius for E = 1. 

Eg (column 6): value of E„ required for R mags: to equal the observed radius, 
rfmax (last column): maximum distance for a single SN explosion with E = 1. 



Table 4. Constraints from radio continuum+ 



Stellar type 


K.41 


S 1420 


ft 






(s- 1 ) 


(mJy) 


(pc) 


(pc) 


03V 


741 


3580 


141 


198 


06V 


219 


1060 


42 


173 


B0V 


14.5 


70 


3 


161 



For a distance of 15 kpc. 



integrate from to R m : dS u (r) = - Jj^'" 4 n r 2 n 2 e a (1) dr, 
where S u (r) is the rate of ionizing photons at distance r, N u 
is the total number of ionizing photons emitted by the star, n e 
the electron number density (assumed equal to n s ) and a <2) = 
3.09 x 10~ 13 cm 3 s" 1 is the recombination coefficient excluding 
captures to the ground level. Solving for R m , we obtain 



R m =R i (l+<R/R i ) U \ <R = 



3N U 



Annla^R 2 ' 



(1) 



If R m > R„, then the HI shell is fully ionized (and should 
not exist!) and some photons escape freely. For our parameters, 
H = 0. 1 9 Af Mi 47 pc, where N Ut 4j is the ionizing photon rate in units 
of 10 47 s -1 . The total number N u of UV ionizing photons also 
fixes the total observed flux density S v at a given frequency and 
is given by (e.g., ??) 



N„ = 0.76 x 10 47 r 4 45 vg^Sv d 2 



kpc* 



(2) 



where T4 is the electron temperature in units of 10 4 K, vghz = 
1.42 the frequency in GHz and S v is in Jy. Using T4 = 0.7 and 
dkpc = 15, we obtain S 1.42 = 4.8 N u ,4i mJy. This estimate for 
S v is an upper limit, because photons can escape through inho- 
mogeneities of the shell (on scales smaller than the beam) or 
through the opening in the north. 

Table [4] shows the results of these simple calculations for 
three different stellar types, chosen to illustrate different possi- 
bilities. Evidently the observations are totally incompatible with 
the presence of an 03V star (or any star with the same value of 
N„) as the HI shell would be essentially totally ionized. In order 
for the neutral shell to be actually present, there would then have 
to be a substantial amount of ionized gas within Ri to absorb 
the ionizing photons. This goes against the canonical view of a 
shell or supershell as surrounding an essentially empty cavity. 
Figure[3]also fails to show any significant amount of continuum 
emission in excess of the diffuse emission seen in this general 
direction which, given the large distance of the shell, is likely to 
be foreground emission. 



At the other extreme, a B0V star (or any star with N u < 
1 .5 x 10 48 s~' would ionize a negligible part of the H I shell, here 
hardly 1 pc. We can use eq. ([D to estimate the ionizing photon 
rate, which would be enough to ionize a thin layer 1 beam size 
in thickness (1 ' or about 5 pc). This gives logYV„ = 48.90, cor- 
responding to an 08.5V star. Note that WR stars have log N u in 
the range 48.6 to 49.4 (Crowther 2007). 

As for the intermediate case of an O6V star, Table [4] shows 
that it would ionize about 13 pc of the inner HI shell. This 
ionized gas would have an angular thickness of about 3 ' and 
a flux density of about 1 Jy. An upper limit to the brightness 
temperature Tb can be obtained by assuming this flux to origi- 
nate from an annulus 35 ' in radius and 3 ' in thickness (giving 
an area of about 660 arcmin 2 or 390 CGPS beam areas at 1420 
MHz). We obtain T B < 2.3 mJy /beam w 0.4 K or, given the 
measured noise of 0.082 K, about 5 cr. An alternative means of 
estimating Tb is to use an approximate emission measure given 
by EM = j n 2 e dl ~n 2 e R ~ 100 pc cm , from which the bright- 
ness temperature is Tb = (EM/566 pc crrT 6 ) K = 0.17 K. Given 
the spatially varying continuum surface brightness within the in- 
ner shell boundary, such a faint ionized layer would not be easily 
detected by our observations. 

Summing up, if the HI shell is relatively homogeneous, 
which is far from certain, we can rule out the presence of 
any star more luminous than O6V or with more than about 
2 x 10 49 ionizing photons s _1 . Placing the shell at the slightly 
closer distance of 1 1 .5 kpc would have the effect of making it 
slightly thinner and closer to any central star with the result that 
for any given spectral type, more of the shell would be ionized. 

4.3.2. A clumpy shell 

In addition to an obvious opening in the north, the H I shell does 
appear to be inhomogeneous on scales comparable to the beam 
size (about 1 ' or 4.5 pc at a distance of 15 kpc). If it were inho- 
mogeneous on scales smaller than the beam (thus undetectable 
in the CGPS image), the H I shell could let a significant num- 
ber of photons escape, thus leading to a more diluted and fainter 
ionized layer. 

The existence of inhomogeneities, whether they are called 
clumps, filaments, clouds, or cloudlets, has been invoked in a 
number of contexts. As a possibly extreme case (with regards to 
size), ? recently postulated the existence, within the stellar wind 
cavity of the 2'.2 nebula N49, of a few hundred cloudlets of inter- 
nal density 10 5 crrT 3 and of a radius as small as 0.05 pc. These 
would provide the dust needed to explain the 24 yum emission ob- 
served within the N49 cavity through erosion and evaporation. 
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We have no way of knowing the degree or type of inho- 
mogeneities which could be present, however, for illustrative 
purposes, we consider inhomogeneities in the form of clumps 
with a radius r — 1 pc and internal density n' = 10 2 cm~ 3 , 
thus each clump would have a mass m = 10.4 m Q . With a 
total HI shell mass of 2.3 x 10 5 m Q , this implies that there 
are N = 2.2 x 10 4 such clumps. The volume filling factor is 
/ » {AnPN l "i)l{4nR 2 i AR) = r 3 N/(3R 2 AR) ~ 7.6 x 1(T 3 . To es- 
timate the transparency of the shell, we note that as viewed from 
the center of the shell, the clumps cover a fraction e of the shell 
area where e ~ Nnr 2 /(4jtR 2 ) = 3fAR/(4r) * 0.23. 

In other words, about 77% of the ionizing photons would es- 
cape through the shell or, taking the case of an 03V star, only 
168 x 10 47 ionising photons s _1 could contribute to partially ion- 
ize the clumps from Ri to R„, resulting in an undetectable radio 
continuum surface brightness. 

In summary, if the H I shell is basically homogeneous, there 
can be no star producing more than about 2 x 10 49 ionizing 
photons per second. Depending on the degree of dumpiness of 
the shell however, such powerful stars could still be present. 
Nevertheless, given the possibly large kinematical age of the 
shell (up to 17 Myr), it is most likely that there are no O stars 
left inside the shell and that any B star, if present, contributes 
negligibly to ionization. 

4.4. Triggered star-formation? 

Shocks in expanding supershells are widely believed to be the 
primary mechanism for triggering star formation (?). Shells 
behind shock fronts experience gravitational instabilities that 
may lead to the formation of large condensations inside the 
swept-up material, and some of these may produce new stars 
(?). An increasing body of observational evidence confirms the 
importance of this mechanism (e.g. ????). Given the size of 
GSH91.5+2-114, it seems reasonable to expect some signs of 
recent star-formation activity in or near the shell border. 

4.4.1. Radio continuum sources 

The 1420 MHz radio continuum image (Fig. [3} shows that a few 
extended, yet relatively compact objects are found near the in- 
ner periphery of the H I shell. Apart from G9 1.11 + 1.57, the so- 
called "head" of the HII complex BG 2107+49 discussed by 
?, there are only two such objects within the boundary of the 
H I shell. Could these be H II regions whose formation was trig- 
gered by the expanding H I shell? An unambiguous answer to 
this question requires a determination of both the distance and 
the radio spectral index. 

We can obtain an estimate of the spectral index a (S v oc v~ a ) 
by measuring the flux densities at 1420 and 408 MHz. Yet for 
both G9 1.56+0.97 and G92.24+1.57, a substantially bright com- 
pact radio source is present in their immediate vicinity (named 
S i and 5 2, see insets of Fig. [3]). Whereas this poses no problem 
at 1420 MHz, the larger beam at 408 MHz results in a partial 
blend of these point sources with the nearby extended structure, 
particularly severe for G92.24+1.57. 

To obtain the spectral indices, we first subtracted all point- 
like sources from the images at both frequencies with a two- 
dimensional Gaussian-fitting routine (Fig.[3]is the result at 1420 
MHz). Because both S \ and S 2 were successfully subtracted at 
this frequency, we then used the imview prograrrQ to measure 



4 This program is part of the DRAO export software package 



the flux density of the remaining extended source. The error was 
estimated by using slightly different background levels. 

On the 408 MHz image however, the larger beam size (« 
4' x3') results in the Gaussian-fitting routine which finds both S 1 
and S 2 to be extended, in contradiction with the fact that both are 
unresolved even at 1420 MHz. Inspection of the 408 MHz image 
showed however that the program did successfully subtract the 
combined flux density of the extended structure plus the nearby 
compact source, for both G9 1.56+0.97 and G92.24+1.57. In or- 
der to obtain the separate flux densities, we proceeded in two 
stages. Firstly, from the point source flux density at 1420 MHz, 
we estimated the 408 MHz flux density, starting with a trial spec- 
tral index a = 0.75, which is representative of the spectral index 
of extragalactic radio sources between 178 and 1400 MHz (?). 
We then created a Gaussian source with this flux density with the 
408 MHz beam size that we removed from the original image. 
We varied a until the best artefact-free subtraction was found. 
The flux density of the extended structure was then estimated 
as the difference between the combined flux density and that of 
the nearby point source. The 408 MHz flux densities of the two 
point sources, S 1 and S 2, were cross-checked by comparing our 
408 MHz values with the 327 MHz values from the Westerbork 
survey (?). The values for S 1 and S 2 at 327 MHz are 57 mJy and 
90 mJy, respectively, in satisfactory agreement with our determi- 
nations. 

Table [5] summarizes the obtained flux densities and spectral 
indices, together with some independent measurements. Both 
sources S 1 and S 2 have a non-thermal spectrum consistent with 
their origin as extragalactic radio sources. As for G91.56+0.97 
and G92.24+1.57, despite a few discrepant measurements, prob- 
ably arising from a difference in background removal and/or in- 
clusion of nearby point sources, they show a spectral index con- 
sistent with a thermal nature. 

Without a distance estimate it is impossible to ascertain 
whether they are associated or not with the H I shell, however 
given their position on or near the H I shell, the possibility that 
they might be H II regions whose formation was triggered by the 
expanding H I shell is worth pursuing. 

Both sources are unfortunately too faint for a significant ab- 
sorption spectrum to be obtained, which would enable us to set 
a limit on the distance. Another way to try to estimate their 
distances is to look for signatures in the HI emission distri- 
bution around these sources (since unfortunately no molecular 
data is available for this region). An inspection of the entire HI 
data cube shows a minimum in the velocity range from about 
-46.0 to -50.0 km s 1 in the area of G91.56+0.97. A well de- 
fined arc-shape structure of enhanced emissivity surrounds the 
minimum toward lower galactic longitudes. In Fig. [5] the con- 
tour delineating the HII region G9 1.56+0. 97 is shown superim- 
posed on the H I emission distribution averaged in the velocity 
interval mentioned above. The excellent morphological correla- 
tion observed between both structures would put G9 1.56+0. 97 
at a distance of about 6.5 kpc, implying that it is not related 
to GSH91.5+2-114. As for G92.24+1.57 we did not find any 
clear H I structure that could be associated with this source. 

As mentioned previously, the source G91.ll + 1.57, associ- 
ated with the H II complex BG 2107+49 (van der Werf & Higgs 
1990), also lies near the periphery of the large HI shell. Could 
the two be associated? ? obtained an H103or recombination line 
spectrum, which showed this source at a velocity of -79.5 ±1.1 
km s _1 . The resolution of their NRAO 43-m observations was 5'. 
Van der Werf & Higgs (1990) obtained a similar value of-80 + 1 
kms -1 , using DRAO HI absorption spectra and an HI 12a re- 
combination line spectrum obtained with the 100-m Effelsberg 
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Fig. 5. HI emission distribution averaged between -46.0 and - 
49.3 km s '.The black contour corresponds to the 8.5 K level of 
the 1420 MHz radio continuum emission. 



telescope at a resolution of 2.'6. These observations would seem 
to place G9 1.1 1 + 1.57 and our HI shell at significantly differ- 
ent distances. However, the VLA 5" radio continuum image of 
Higgs et al. (1987) shows the head of BG 2107+49 to consist of 
a set of discrete knots of emission about 10" in size and a more 
extended and diffuse region about 3' in diameter. Although we 
cannot completely rule out the possibility that some of the knots 
could be at the larger distance of our H I shell, the fact that the 
H103a- recombination line spectrum of Higgs et al. (1987) is 
extremely well fitted by a single line at -79.5 km s _1 with es- 
sentially no residuals near - 1 14 km s" 1 would seem to rule this 
out. 

4.4.2. Infrared sources 

We feel that the case for triggered formation is stronger for 
IRAS 21 147+5016. This IRAS point source is seen projected 
inside the shell wall and appears itself surrounded by a smaller 
HI bubble (Fig.Q]). This smaller HI cavity resembles the ones 
found by ? and ?. In both cases the authors found a small cavity 
immersed inside a larger HI shell. However, their interpretation 
about the origin of the structures differed. ? concluded that both 
the large HI shell and the smaller bubble had been created by the 
same star, HD 197406, during different stages of its evolution. 
On the other hand, ? found an HI cavity around the OB associ- 
ation Bochum 7 and concluded that this association might have 
been born as a consequence of the evolution of the large shell 
GS263-02+45. 

From Fig. 1 we can infer that the size of the smaller cav- 
ity is about 0°3. At the distance of GSH91. 5+2-114, this im- 
plies a linear size of about 80 pc. Given that this cavity can 
be observed over at least 10 kms -1 (see Fig. 2), a lower limit 
for the expansion velocity can be assumed as 5 kms -1 . Under 
these conditions we derived an upper limit of about 5 Myr for 
the dynamical age if a stellar wind origin is considered. This 
age is significantly smaller than the one obtained for the large 
shell (7-17 Myr), but larger than the duration of the WC phase 
(10 5 yr), suggesting that the O-phase of the current WC star 



would also have contributed to its formation. Taking into ac- 
count that progenitors with a mass of some 40 to 50 M G are 
suggested for late WC stars (?), the dynamical age estimated 
for the small cavity is consistent with the time that these stars 
stay on the main sequence, between 3.7 and 4.9 Myr (?). Based 
on the age estimate difference, the most probable scenario is one 
where the large shell, GSH 91.5+2-1 14, was created by the joint 
action of several massive stars and where at least one of them 
has already exploded as a supernova. Furthermore, as the shell 
evolved, new stars may have been triggered in its dense border, 
IRAS 21 147+5016 being one of them. It is of interest that this 
smaller bubble (Fig. [3} appears to be opened on the exterior side 
of the larger H I shell. This would be expected for a star whose 
formation would have taken place toward the exterior of the ex- 
panding shell, leading it to first burst outwards. 

5. Conclusions 

The measured kinetic energy of expansion of GSH91. 5+2-1 14 
is far too large to have been produced solely by the action of one 
or all of the four DWCL sources possibly associated with the 
H I shell. That the shell is very symmetrical whereas the DWCL 
sources are located well off-center suggests that if they are in- 
deed inside the shell, their shaping influence is minimal. We have 
found no evidence for the presence of other massive stars, but 
absorption would preclude the detection of these objects at the 
inferred distance. 

An interpretation as the HI shell of a single SNR is 
also not tenable. As for GSH 138-01-94 and GSH 90+03-99, 
GSH 91.5+2-114 has likely been caused by the combined winds 
and SN explosions of a number of massive stars. 

The H I shell appears open to the north, i.e. in a direction 
away from the Galactic plane, suggesting that the shell has burst 
out of the Galactic disk. A filament displaced from the north- 
ern boundary (Filament B, Fig. [TJ could be the remains of the 
top (now missing) part of the shell. The faint radio continuum 
emission, consisting of filaments more or less aligned in a di- 
rection perpendicular to the Galactic plane, lends support to the 
breakout hypothesis. 

Two relatively compact thermal sources (Table |5), seen in 
projection near or on the boundary of the HI shell, could 
have formed in gas compressed by the expanding shell. One of 
the DWCL sources, IRAS 21 147+5016, seen projected inside a 
smaller H I shell, could also be a product of triggered star for- 
mation. We cannot rule out the possibility that the other DWCL 
sources could also have formed in the same manner. 
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Table 5. Flux densities and spectral indices of compact sources + (S y oc v~ a ) 



# 


Object 


^1420 (mjy) 


cr 1420 (mJy) 


F 4 o 8 (mJy) 


o- 4 o 8 (mJy) 


a 


Aa 


Notes" 


1 


S 2 


37 


2 


76 


9 


0.58 


0.10 


This paper, see text 


2 


G91.56+0.97 


191 


19 


159 


27 


-0.15 


0.16 


This paper, see text 


3 




261 


11 


346 


28 


0.09 


0.06 


vdWH, HvdW" 


4 




238 


9.5 


79.2 


1.1 


-0.88 


0.03 


Kerton et al. (2007) 


5 


Si 


25 


2 


59 


4 


0.69 


0.08 


This paper, see text 


6 


G92.24+1.57 


132 


10 


104 


15 


-0.19 


0.13 


This paper, see text 


7 




99.4 


6.9 


168 


31 


0.42 


0.16 


Kerton et al. (2007) c 


8 




146 


22 










vdWH 



+ Sources 5 1 and 5 2 are the nearby point sources seen in the insets of Fig. 3 

" HvdW = Higgs & van der Werf (1991), vdWH = van der Werf & Higgs (1990) 

b Flux density at 408 MHz includes contribution from source S 2 (29P62 at 1420 MHz - HvdW). Note that the spectral index listed there 

from a linear fit using measurements at five frequencies between 408 MHz and 4.85 GHz 
c Likely includes source 5 1 at 408 MHz 



